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ABSTRACT 



^ ' Context. We present the observational results of multi-colour optical monitoring of eight red blazars from 2003 September to 2004 February. 
Aims. The aim of our monitoring is to investigate the spectral variability as well as the flux variations at short and long time scales. 
Methods. The observations were carried out using the 1.0 m robotic telescope of Mt. Lemmon Optical Astronomy Observatory, in Arizona, 
USA, the 0.6 m telescope of Sobaeksan Optical Astronomy Observatory and the 1.8 m telescope of Bohyunsan Optical Astronomy Observatory, 
in the Republic of Korea. 

Results. During the observations, all sources show strong flux variations with amplitudes of larger than 0.5 mag. Variations with amplitudes of 
over 1 mag are found in four sources. Intraday variations with amplitudes larger than 0.15 mag, and a rapid brightness increase with a rate of 
~ 0.2 mag per day in four days, are detected in S5 0716+71. We investigate the relationship between the colour index and source brightness for 
Q ' each source. We find that two out of three FSRQs tend to be redder when they are brighter, and, conversely, all BL Lac objects tend to be bluer. 

, In particular, we find a significant anti-correlation between the V- 1 colour index and R magnitude for 3C 454.3. This implies that the spectrum 
Q I became steeper when the source was brighter, which is opposite to the common trend for blazars. In contrast, significant positive correlations 
5-H are found in 3C 66A, S5 0716+71, and BL Lac. However, there are only very weak correlations for PKS 0735+17 and OJ 287. 
C/j Conclusions. We propose that the different relative contributions of the thermal versus non-thermal radiation to the optical emission may be 
responsible for the different trends of the colour index with brightness in FSRQs and BL Lac objects. 
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■ 1. Introduction tribution (SED), with two broad spectral components. The first, 

lower frequency component is generally interpreted as being 

Blazars, including BL Lac objects and flat-spectrum radio due to synchrotron emission, and the second, higher frequency 

quasars (FSRQs), are the most extreme class of active galac- ^^-^^ -^^^^^^ Compton emission. Although dif- 

tic nuclei (AGNs), characterized by strong and rapid variabil- f^^g^^ ^^^^^^^ j^^^g diff^erent peak frequencies, the two peaks in 

ity, high polarization, and apparent superluminal motion. These ^^^^-^ ^^^^ separated by approximately flie same amount, 

extreme properties are generally interpreted as a consequence j ^ g.^Q ^^^^^^^ frequency (Fossati et al. 1998). According 

of non-thermal emission from a relativistic jet oriented close ^^ffg^g^j p^^,^ frequencies, blazars are divided into two 

to the line of sight. As such, they represent a fortuitous natural subclasses: the low-energy-peaked blazars (red blazars), which 

laboratory with which to study the physical properties of jets, j^^^g synchrotron peaks in the IR/optical range, and the high- 

and, ultimately, the mechanisms of energy extraction from the energy-peaked blazars (blue blazars), which have synchrotron 

central supermassive black holes. p^^j^^ UV/X-ray energies. All of the TeV y-ray-loud AGNs 

In the field of AGN observations, one of the most impor- are blue blazars, while aU strong GeV y-ray AGNs are red 

tant discoveries in the last decade has been that blazars emit blazars (Mattox et al. 1997; Hartman et al. 1999; Bai & Lee 

a substantial fraction, sometimes most of their power, at y-ray 2001). 
energies (GeV and TeV, e.g. Ulrich, Maraschi & Urry 1997; 

Catanese & Weeks 1999). The y-ray emission of blazars indi- Although a point of agreement about the second compo- 

cates a double-peak structure in the overall spectral energy dis- nent is that the y-rays of blazars are produced in relativistic jets 

by inverse Compton scattering, the origin of the seed photons 

Send ojfprint requests to: M. F. Gu, e-mail: gumf@shao.ac.cn (optical/IR), the location and size of the emitting region, and 

* Table 2 is only available in electronic at the CDS via the degree of relativistic beaming of the high-energy radiation, 

anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via are all stifl unknown. Seed photons may come from the jet it- 

,http://cdsweb.u-strasbg.fr/Abstract.html_ self (Synchrotron Self-Compton model), from an accretion disk 
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around a supermassive black hole at the base of the jet, or else 
from photons of the broad emission line region (e.g. Wehrle et 
al. 1998; CoUmar et al. 2000). Some models suggest that the 
flares of y-rays in blazars are caused by the increase and de- 
crease of the bulk Lorentz factor, while oilier models suggest 
changes in both the injection rate and the spectral shape of the 
injected electrons (Mukherjee et al. 1997; Hartman et al. 2001; 
Plan et al. 2002; ViUata et al. 2002). 

Simultaneous multiwavelength observations have shown 
that both TeV and GeV y-ray blazars vary much more vio- 
lently at and above their peaks, than below these peaks, and that 
flux variations at these two peaks are correlated, indicating that 
both components are produced by the same electron population 
in the relativistic jet (e.g. Ulrich et al. 1997; Sambruna et al. 
2000; Hartman et al. 2001). Since red blazars have synchrotron 
emission peaking in the IR/optical bands, observations at these 
wavelengths are thus crucial in understanding the nature of y- 
ray emission in red blazars. The near-lR/optical bands are ideal 
for monitoring the flux and spectral variabihty, as well as to de- 
termine the frequencies and frequency shifts of the synchrotron 
peaks for red blazars. The flux and spectral variability, together 
with the synchrotron peak frequency and its shift, aUow one to 
derive information on the radiating relativistic electrons, and 
on the emission region (Ghisellini et al. 1997). Moreover, the 
difl'erences in the variabihty behaviour, in the IR/optical bands, 
of non-y-ray-detected red blazars, as compared to those of the 
y-ray-loud ones, may strongly constrain models for y-ray emis- 
sion. 

The study of the properties of the SED is an excellent diag- 
nostic tool for theoretical models. However, the availability of 
simultaneous multiwavelength observations is generaUy inade- 
quate in tracing the time evolution for both the synchrotron and 
inverse Compton components. A relatively simple method is to 
investigate time evolution in a selected region of the spectra, 
such as in the optical. Although the optical region is narrow 
with respect to other spectral regions, it may provide a rela- 
tively large amount of information. In this case, the possible 
presence of other components in addition to the synchrotron 
emission, e.g. thermal emission from the accretion disk around 
the central engine, may complicate the situation. Nevertheless, 
statistical analysis of the optical spectral variability could en- 
able us to constrain theoretical models; as claimed by Vagnetti, 
Trevese & Nesci (2003), the spectral variability, even when re- 
stricted to the optical band, can be used to set limits on the rela- 
tive contribution of the synchrotron component and the thermal 
component to the overall SED. 

In this paper, we investigate the spectral variability, as well 
as the multi-band flux variability for eight typical red blazars, 
which consist of three FSRQs and five BL Lac objects. The 
spectral variability of the FSRQs is compared with that of the 
BL Lac objects. Based on analysis of the spectral variability, 
we attempt to derive information to constrain the theoretical 
models. The observations and data reduction are given in Sect. 
2. The flux variability is shown in Sect. 3, and the spectral vari- 
ability is described in Sect. 4. The discussions are given in Sect. 
5, and the main results are summarized in Sect. 6. 



2. Observations and data reduction 

Multi-site photometric observations of eight blazars were car- 
ried out for a combined total of 50 nights, from 2003 September 
to 2004 February, using the 1.0 m robotic telescope of Mt. 
Lemmon Optical Astronomy Observatory (denoted as L in 
Table 1; Han et al. 2005), in Arizona, USA, for 33 nights, the 
0.6 m telescope of Sobaeksan Optical Astronomy Observatory 
(denoted as S in Table 1) and the 1.8 m telescope of Bohyunsan 
Optical Astronomy Observatory (denoted as B in Table 1), in 
the Republic of Korea, for 14 and 3 nights, respectively. The 
observations were made in 5 sessions to investigate both long- 
term and intraday variations: 5 nights in September, 10 nights 
in October, 13 nights in November, 16 nights in December, 
2003, and 6 nights in February, 2004. All telescopes were 
equipped with CCD cameras, and UBV (Johnson) and RI 
(Cousins) filter sets. LN2 cryogenic systems were used for both 
the 1.8 m and 0.6 m telescopes, while a thermoelectric cooling 
system was used for the 1 .0 m telescope. These three telescopes 
cover different field of views: 22.2x22.2 arcmin, 20.5x20.5 ar- 
cmin, and 11.7x11.7 arcmin for the 1.0 m, 0.6 m, and 1.8 m 
telescopes, respectively. Considering the aperture sizes of the 
telescopes and the seeing conditions, difl'erent exposure times 
were applied. While observing, we tried to position each target 
object at the same location on the CCD surface, to within a few 
pixels, in order to achieve more eflicient and accurate photom- 
etry. Sky flat images were taken at both dusk and dawn when 
available. Bias and dark images were taken after obtaining sky 
flat images. 

All images were preprocessed with the standard proce- 
dures, using the IRAP' software, including bias and dark sub- 
traction, and flat-fielding. For images obtained from the 0.6 m 
telescope, the fringe pattern was carefully corrected. We used 
the APPHOT package in IRAF to measure the instrumental 
magnitudes of the blazars and comparison stars. The seeing 
variations during one night may cause variations of the in- 
strumental magnitude, and there is a trend of the instrumental 
magnitude variations with the FWHM variations of a star from 
night to night (Clements & Carini 2001). In order to correct the 
seeing variation effect during one night, and to try to include the 
'total' flux from objects when seeing varies (especially when 
seeing is poor), we set the radius in the aperture photometry to 
be proportional to the FWHM of bright, well-exposed stars in 
the image frames containing the objects, during one night (Lee 
et al. 2003). After experimenting with various aperture sizes, 
we set the aperture radius at 4xFWHM, the inner radius of the 
sky annulus at 5xFWHM, and the width of the sky annulus at 
10 pixels, to maximize the S/N ratio. 

In Table 1, the source lAU name and the alias name are 
given in Cols. 1 and 2, respectively. The classification, fil- 
ter sets, telescopes, and the detection in the y-region for each 
source are shown in Cols. 3-6. The source magnitude is cali- 
brated with respect to the standard stars, for which references 
are given in Col. 7 of Table 1. The observational errors are es- 

' IRAF is distributed by the National Optical Astronomy 

Observatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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timated from the rms differential magnitude between the cali- 
bration star and another standard star used for checking, 



where m, = {mc - tnK)i is the differential magnitude of the 
calibration star C and check star K, while m = mc - mk is the 
differential magnitude averaged over the entire data set, and 
is the number of observations on a given night. For each object, 
the typical rms error is between 0.01 and 0.02 mag. 

3. Flux variability 

The results of our monitoring program are listed in Table 2, 
available in electronic form at the CDS. Column 1 is the object 
name. Col. 2 is the Juhan Date, Col. 3 is the magnitude. Col. 4 
is the rms magnitude error. Col. 5 is the filter used, and Col. 6 is 
the observatory. For each object, the historical photometric re- 
sults are briefly presented. The overall magnitude variations of 
each object are investigated. Variations on time scales of days, 
and on shorter time scales (e.g. hours), in several objects are 
also reported. 

3.1. 3C66A 

The blazar 3C 66A is a member of the Einstein Slew Survey 
sample (Elvis et al. 1992). Both GeV (Dingus et al. 1996; 
Mukherjee et al. 1997) and TeV (Neshpor et al. 1998) y- 
rays have been detected in this highly polarized (Mead et al. 
1990) blazar. This source has displayed infrared, optical, ul- 
traviolet, and X-ray variability on different time scales (Ghosh 
& Soundararajaperumal 1995; De Diego et al. 1997; Carini, 
Noble, & Miller 1998). 3C 66A was the target of an extensive 
multiwavelength Whole Earth Blazar Telescope (WEBT) mon- 
itoring campaign from July 2003 through April 2004, involving 
observations in radio, infrared, optical. X-ray, and very-high- 
energy y-ray bands (Bottcher et al. 2005). At all optical fre- 
quencies, a gradual brightening of the source over the course of 
the campaign was observed, with a maximum at ^ » 13.4 on 
Feb. 18, 2004. Microvariability with flux changes of ~ 5 % on 
time scales as short as ~ 2 hr, and several bright flares on time 
scales of several days, were detected (Bottcher et al. 2005). 
Our observational data on this source have been included in 
the WEBT campaign paper (Bottcher et al. 2005). 

Our BVRI light curves are shown in Fig. 1. Data obtained 
from different observatories are represented by different sym- 
bols. During our observing runs, the overall magnitude varia- 
tions were ^B = 0.59, AV = 0.55, Mi = 0.50, and A/ = 0.52 
over 88 days. Our observations of this object consist of four 
sessions: September, October, November, and December 2003. 
From Fig. 1, it can be seen that the source stays at a relatively 
low state, R ~ 14.3 in the September session. However, 3C 66A 
brightened io R ~ 14.0 in the following three sessions. It can 
also be seen that the source brightens gradually in the October, 
and December sessions, but fades in the November session, al- 
though there are fluctuations in each observing run. We find 
that the brightest magnitude of October, R ~ 13.8 is same as 
that of December. The time separation of these two epochs is 



about 57 days. Short- time- scale variations (e.g. intraday varia- 
tions) were also searched for during our observing runs, but no 
significant such variations were found. 

3.2. PKS 0420-01 

This quasar has revealed a strong variability in the optical band 
(Webb et al. 1988). A noticeable flare was detected in late 1979, 
when a 1 .3 mag increase in 5 days was registered, followed by a 
1.7 mag decrease in 23 days (Villata et al. 1997). A strong vari- 
ability is also evident from the 20 yr light curve (from 1970 to 
1990) reported by Smith et al. (1993). Wagner et al. (1995) no- 
ticed flux variations with time scales of the order of 1-10 days. 
The flare of 1992 February-March was the highest optical state 
observed until then (R = 14.6); in that period EGRET regis- 
tered the highest y-ray flux density. Since low fluxes or non- 
detections at y-ray energies correspond to low optical states, a 
direct correlation between the optical and y-ray emission was 
suggested. The most noticeable variation was a fall of 2.64 
mag in 40 days observed from 1995 September 15 to October 
25 in the R band (Villata et al. 1997). Microvariability with 
A/? = 0.12 mag in 40 minutes was also detected (Villata et al. 
1997). 

The VRI light curves of PKS 0420-01 are shown in Fig. 
2. During our observing runs, the overall magnitude variations 
were AV = 1.41, AR = 1.48, and Al = 1.54. The maximum 
R = 15.06 happened on October 30, 2003 (JD=2452942.87), 
and the minimum R = 16.54 was on December 25, 2003 
(JD=2452998.76). During the September run, tiie fight curve 
over three days (September 29, 30, and October 1) shows that 
the source was near R ~ 15.4, and it slightly brightened over 
this period, but with only about 0.12 mag variation in this band. 
During the October session, the source gradually brightened 
from R = 15.74 + 0.03 on October 26, to R = 15.06 ± 0.01 
on October 30, which is the maximum over the whole mon- 
itoring period. The average brightness during these five days 
is R ~ 15.4, which is same as the average magnitude of the 
September session. We find that the source also stays at the 
same brightness level (R ~ 15.4) in the November run, al- 
though there are fluctuations from day to day. The R magnitude 
was ~15.47 on November 20, and it brightened to ~ 15.25 on 
November 21 . On November 22, the source faded to R ~ 15.43, 
and after that, it stayed atR ~ 15.4. However, the source had 
faded by about 1 mag toR ~ 16.5 in the December session. 

3.3. S5 0716+71 

The BL Lac object S5 0716-1-71 was included in the S5 cata- 
log of the Strong Source Survey performed at 4.9 GHz (Kiihr 
et al. 1981). The largest optical database on S5 0716-1-71 ever 
pubUshed was collected and analyzed in Raiteri et al. (2003), 
which consists of a total of 4854 UBVRI data points obtained 
from eight observatories from 1994 to 2001. Four major opti- 
cal outbursts were observed: at the begiiming of 1995, in late 
1997, at the end of 2000, and in autumn 2001. In particular, an 
exceptional brightening of 2.3 mag in 9 days was detected in 
the R band just before the BeppoSAX pointing of October 30, 
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Table 1. List of the eight blazar sources in the optical monitoring sample. 



Object 


Name 


Class. 


Filters 


Observatory 


y-detected 


Comparison Stars 


0219+428 


3C66A 


BL Lac 


BVRI 


L,S 


Yes 


Fiorucci & Tosti (1996)" 


0420-014 


PKS 0420-01 


FSRQ 


VRI 


L,S 


Yes 


Smith feBalonek (1998) 


0716+714 


S5 0716+71 


BLLac 


BVRI 


B,L,S 


Yes 


Villata etal. (19983)" 


0735+178 


PKS 0735+17 


BL Lac 


VRI 


L,S 


Yes 


Smith etal. (1985) 


0851+202 


OJ287 


BL Lac 


VRI 


L,S 


Yes 


Fiorucci & Tosti (1996) 


1641+399 


3C 345 


FSRQ 


VRI 


L,S 


No 


Smith etal. (1985)" 


2200+420 


BL Lac 


BLLac 


VRI 


L,S 


Yes 


Fiorucci & Tosti (1996)" 


2251+158 


3C 454.3 


FSRQ 


VRI 


L,S 


Yes 


Fiorucci et al. (1998) 



" /-data from Ghisellini et al. (1997) 

* Magnitudes from Gonzalez-Perez et al. (2001) 



Table 2. The observational results. 



Name JD Magnitude Error Filters Observatory 



14.6 - 




JD - 2452900 

Fig. 1. Light curves of 3C 66A in the B,V, R, I bands. The circles represent the data from the LO m Mt. Lemmon telescope, and 
the triangles are from the 0.6 m Sobaeksan telescope. 



2000. They found that the long-term trend shown by the opti- 
cal light curves seems to vary with a characteristic time scale of 
about 3.3 years. The source has been recently monitored at ra- 
dio and optical wavelengths by more than 40 telescopes in the 
northern hemisphere, during a WEBT campaign lasting from 
September 2003 to June 2004. Nesci et al. (2002) presented a 
study of the optical intraday variability (IDV) of S5 0716+71 
over 52 nights. They found typical variation rates of 0.02 mag 
per hour, and a maximum rising rate of 0.16 mag per hour. 
S5 0716+71 has been observed with INTEGRAL on 2-6 April 
2004 (Pian et al. 2005). It was detected with IBIS/ISGRI up to 



60 keV, with a flux of ~ 3 x 10"" erg s'^ cm-^ in the 30-60 
keV interval. 

S5 0716+71 was the most extensively monitored source 
in our observations. The observations span about 140 days, 
occurring in September, October, November, December 2003, 
and February 2004. The light curves in the Johnson's BY and 
Cousins' RI bands are shown in Fig. 3, where the different sym- 
bols refer to the various telescopes listed in Table 1. It can be 
seen from Fig. 3 that the source exhibited strong variability, 
with similar trends (but with slightly different amplitudes) in 
all bands. The overall magnitude variations were AB = 1.42, 
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H 1 1 1 1 1 1 H 



50 

JD - 2452900 



Fig. 2. Light curves of PKS 0420-0 1 in the V, R, I bands. The symbols are the same as in Fig. 1 . 



AV = 1.37, AR = 1.37, A7 = 1.36, where the differences 
are mainly due to the different temporal samplings. During our 
observing runs, the maximum /?-band brightness, R = 12.62, 
was recorded on 10 December 2003 (JD=2452983.8756). We 
note that a rapid brightness increase was detected during the 
4 days of 19 -23 November 2003, and which is shown in Fig. 
4. The overall variation in the R band was AR = 1.15, from 
R = 13.94 ± 0.01 on JD=2452962.916 to = 12.79 + 0.01 
on JD=2452966.854, both observed from the Mt. Lemmon 
Observatory. 

Interestingly, we find that the shape of the profile of this 
rapid brightness increase is significantly Unear. We have fit- 
ted least-square regression lines to the light curves in the BVRI 
bands obtained from all three telescopes. The slopes of the best- 
fit lines are given in Col. 2 of Table 3, together with the regres- 
sion coefficients for the best fits in Col. 3. It is noted that the 
slopes are essentially the same for the different passbands. The 
differences in the slopes are always < Icr, and are therefore not 
statistically significant. Our results show that the rapid bright- 
ening can be described by linear trends on a magnitude scale 
(which, therefore, corresponds to exponential intensity varia- 
tions), and this rapid brightening rate was ~ 0.2 mag per day 
over 4 days. However, it should be noted that this rapid bright- 
ening is just a mean behaviour, and may also be due to the lack 
of data in the gaps. 

Rapid variations of the flux of S5 0716-1-71 have been ob- 
served on different occasions (e.g. Ghisellini et al. 1997; Sagar 
et al. 1999; Xie et al. 2004). In our observations, we also 
searched for intraday variations. Variations, with amplitudes 
of >0.15 mag over a few hours in one night, were found on 
November 23 and December 10, 2003, and which are shown 
in Fig. 5, and Fig. 6, respectively. Both of them were detected 



with the Mt. Lemmon Observatory. The lower panels in these 
figures show the difference between the instrumental magni- 
tude of the two field standard stars. In Fig. 5, it can be seen that 
the source brightened from y = 13.25±0.02(/? = 12.89±0.01) 
to y = 13.17 + 0.02 (R = 12.79 + 0.01) within 48 min, and af- 
terwards it faded to y = 13.45 ± 0.02 (R = 13.06 ± 0.01) in the 
remaining observing time (~ 4.6 hr) for this night. The over- 
all magnitude variations were Ay = 0.28 mag (AR - 0.27 
mag) during this night. The variability parameters C (intro- 
duced by Romero, Cellone & Combi 1999; if C>2.576, then 
the source is variable) are 3.64 (V band), and 5.27 (R band). 
The lower panels show that the maximum deviation of the two 
field standard stars is Ay^ax = 0.07 mag in the V band, and 
ARmax = 0.06 in the R band. In Fig. 6, the source brightened 
fromy= 13.18 + 0.01(5= 13.64 ± 0.01) to y = 13.02 + 0.01 
(B = 13.47 ± 0.01) within about 1.4 hr, then faded. The ris- 
ing rate, of 0.002 mag per minute in both B and V bands, was 
found over the duration of the first 1 .4 hr. This is in good agree- 
ment with the results of Villata et al. (2000), where there was 
a monotonic brightness increase from B = 14.13 ± 0.01 to 
B - 13.98 ± 0.01 in about 130 minutes, with the steepest (hn- 
ear) part having a rising rate of 0.002 mag per minute, and a 
duration of about 45 minutes. Moreover, the authors suggested 
that such a gradient seems to be typical of fast variations, with 
no evident difference between rise and fall. 



3.4. PKS 0735+17 

PKS 0735+17 was first classified as a BL Lac object by 
Carswell et al. (1974). It is both radio (Kiihretal. 1981) and X- 
ray selected (Elvis et al. 1992). Falomo & UMch (2000) gave 
a lower limit on the object's redshift, of z > 0.5, and Xie et 
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Fig.3. Light curves of S5 0716+71 in the B, V, R, I bands. The circles and triangles are the same as in Fig. 1, and the crosses 
are the data from the 1 .8 m Bohyunsan telescope. 

Table 3. Slopes and regression coefficients (r) of the least-square linear fits to the BVRI magnitude variations of S5 0716-1-71, 
during 19-23 November 2003. 



Filters 


slopeicr 


r 


B 


-0.197±0.013 


-0.90 


V 


-0.194±0.013 


-0.90 


R 


-0.198+0.012 


-0.92 


I 


-0.179+0.012 


-0.89 



al. (2002) have estimated its black hole mass. Tommasi et al. 
(2001) observed the optical polarimetry of PKS 0735+17 in 
1999 December, and gathered together previous polarimetric 
data. The source has shown very different levels of polarization 
percentage in past years, from around 1% to more than 30%. 
PKS 0735 + 17 is also known as an optical and infrared intraday 
variable blazar (Massaro et al. 1995; Heidt & Wagner 1996; 
Bai et al. 1998). With observational data spanning over 90 yr, 
the long-term periodicity of this blazar has been investigated; 
its optical variability shows the possibihty of a 4.89 yr period 
(Webb et al. 1988; Smith et al. 1988), and also a 14.2 yr period 
(Fanetal. 1997). 

Our VRI light curves are shown in Fig. 7. During our ob- 
serving runs, the overall magnitude variations were AV - 0.57, 
AR = 0.54, and Al = 0.46. From Fig. 7, it can be seen that this 
source is variable almost all the time, although the overall vari- 
ations are not relatively large. This source was relatively stable 
from September 30 to October 1, 2003. The average R magni- 
tude of PKS 0735+17 was 15.13 + 0.02 on September 30, and 
15.17 ± 0.04 on October 1 . From October 26 to October 30, the 
apparent fluctuations are present in all bands (see Fig. 7). On 



October 30, the maximum R brightness of whole monitoring 
period was reached at 14.81 ± 0.01 (10=2452942.998). In the 
R band, a brightening from 15.26 + 0.01 on November 21 to 
15.09 ± 0.04 on November 23 was detected. One month later, 
the object was observed atR = 15.27 ± 0.02 on December 22, 
and brightened to 7? = 15.16 + 0.02 after two days. FinaUy, it 
faded to = 15.31 ± 0.02 on December 25. 

3.5. OJ287 

OJ 287 is one of the most extensively observed BL Lac ob- 
jects. It has been observed for over 100 years, providing a very 
good historical light curve. In 1972, it reached maximum light, 
with V ~ 12. Pursimo et al. (2000) presented intensive opti- 
cal, infrared, and radio monitoring of OJ 287, taken between 
1993 and 1998. These monitoring results show that the opti- 
cal and infrared light curves displayed continuous variability 
on timescales ranging from tens of minutes to years. The his- 
torical light curve shows several outbursts with a recurrence 
period of about 1 1.65 yr (Sillanpaa et al. 1988a). The 12 yr cy- 
cle was discussed in many occasions (e.g. Villata et al. 1998b; 
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Fig. 4. The variability of S5 0716+71 during 19 - 23 November 2003, in the B, V, R, I bands. The symbols are the same as in 
Fig. 3. 



Abraham 2000; Katajainen et al. 2000; Valtaoja et al. 2000). 
Sillanpaa et al. (1988a) modeled these periodic outbursts by 
using a binary supermassive black hole system, and predicted 
that the next outburst would occur during late 1994. This pre- 
diction was confirmed by several studies (Kidger et al. 1995; 
Sillanpaa et al. 1996a, 1996b; Arimoto et al. 1997; Fan et al. 
1998; Jia et al. 1998; Pursimo et al. 2000; Efimov et al. 2002). 

Our VRI light curves are shown in Fig. 8. During our 
observing runs, the overall magnitude variations were AV = 
1.39, AR = 1.40, and A/ = 1.27. The maximum R = 
14.23 ± 0.01 occurred on September 30 (10=2452913.001), 
whereas the minimum R = 15.63 + 0.02 was on November 23 
(JD=2452966.982). The source was observed at R = 14.23 ± 
0.01 (V = 14.65 ± 0.03, / = 13.55 ± 0.05) on September 30, 
2003, which was recorded at the Mt. Lemmon Observatory. 
About one month later, it faded to R = 15.42 ± 0.02 on 
October 26. After that, it gradually brightened, and reached R = 
14.95 + 0.02 on October 30. In November, the source bright- 
ness stayed at a similar level as that in October. On November 
20, the R magnitude was recorded at 15.43 ± 0.01, but it bright- 
ened to 15.15 ± 0.01 on November 21. Later, the source faded 
to 15.58 ± 0.02 on November 23. In December, a fading from 
R = 15.20 ± 0.01 on December 22, to R = 15.44 ± 0.03 on 
December 25 was detected over three days, at the Mt. Lennmon 
Observatory. 

3.6. 3C 345 

3C 345 was classified as an OV V quasar by Penston & Cannon 
(1970). The typical behaviour for this object has been to show 



about 2 mag outbursts, which occur quite frequently. Such large 
outbursts happened in 1967/1968, 1971/1972, 1982/1983, and 
1991/1992 (Schramm et al. 1993). Sillanpaa et al. (1988b) 
measured V = 14.56 in 1982, while Schramm et al. (1993) 
observed the V-band brightness near 18 mag in 1989. In 
1991/1992, a large 2.5 mag outburst was observed (Schramm 
et al. 1993). Kidger & Takalo (1990) observed the historical 
minimum at B = 18.66. 

The VRI Ught curves are shown in Fig. 9. During our ob- 
serving runs, the overall magnitude variations were AV = 0.66, 
AR = 0.55, and Al = 0.63. The maximum R brightness R = 
15.81+0.02 was observed on December 25 (JD=2452999.055), 
and the minimum R = 16.36 ± 0.01 was on September 30 
(JD=2452912.710). Variations with small ampUtudes on time 
scales of days are found in all bands, in September, October, 
and December 2003. 

3.7. BL Lac 

BL Lac is the prototype of the BL Lac object class of AGNs, 
and one of the best studied blazars. Shen & Usher (1970) have 
investigated its historical light curves, and found a range of 
variation of 4.2 magnitode in the V band, and a strong out- 
burst in which BL Lac varied over nearly its entire range 
(12.4< B <16.7 mag) during 400 days. Investigations of its 
long-term variability have been done by Webb et al. (1988), 
and Carini et al. (1992), which show that the faintest magni- 
tudes were B - 17.99 mag, and V - 16.73 mag. Recently, 
Villata et al. (2004) presented UBVRI light curves obtained 
by the WEBT from 1994 to 2002, including the last, extended 
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Fig. 5. Intraday variations of S5 0716+71 on 23 November 2003, in the V band (left panel) and in the R band (right panel). All 
data are from the 1.0 m Mt. Lemmon telescope. 



BL Lac 2001 campaign. Their analysis of the colour indices 
reveals that the flux variability can be interpreted in terms of 
two components: longer-term variations occurring on a few- 
day time scales appear as mildly-chromatic events, while a 
sttong bluer-when-brighter chromatism characterizes very fast 
(intraday) flares. They suggested that Doppler factor variations 
on a "convex" spectrum could be the mechanism accounting 
for both the long-term variations, and their slight chromatism. 

Our VRI light curves are shown in Fig. 10. During our 
observing runs, the overall magnitude variations were AV = 
1.06, AR = 1.05, and A/ = 0.94. The maximum R-hmd 
brightness was recorded at R = 13.55 + 0.01 on November 
20 (JD=2452963.580), whereas the minimum was detected at 
R = 14.60 ± 0.01 on October 28 (10=2452940.784). BL Lac 
was observed at R ~ 14.15 + 0.01 on September 29. It was 
slightly fainter in late October. On October 26, the R magni- 
tiide was R = 14.37 ± 0.04. It became R = 14.23 ± 0.05 after 
four days, during which small amplitude fluctuations (< 0.2 
mag) were detected from day to day. In the November run, 
BL Lac first brightened to R = 13.55 ± 0.01 on November 20 
(JD~2452963.580), then it faded to R ^ 14.02 + 0.01 in one 
day. On November 22, a brightening of ~ 0.10 mag in the R 
band (0.11 mag in the V band, and 0.08 mag in the / band), 
from/? = 14.03+0.01 at JD=2452965.648 to = 13.93+0.01 
at JD=245 2965. 696, was detected in about 69 minutes. The R 
magnitude was recorded at 13.86±0.01 on December 22, but it 
dropped to 14.23 + 0.01 after one day. After that, it brightened 
slightly to 14.14 + 0.02 on December 24. Finally, it reached 
14.24 ± 0.02 on December 25. 



3.8. 3C 454.3 

3C 454.3 is a GeV y-ray source, and one of the few blazars 
that have been detected by COMPTEL (Blom et al. 1995). The 
light curve of 3C 454.3 in tiie B band from 1966 to 1979 is 
given by Lloyd (1984). Data in the B band from 1971 to 1985 
are reported by Webb et al. (1988), and a flare with a variation 
of 1.28 mag in 63 days was reported in autumn 1979. From 
1986 November to 1987 January, the average B magnitude of 
the source was 16.56, with a variation of 0.92 mag over 44 days 
(Corso et al. 1988). The results of the monitoring program of 
Villata et al. (1997) show that a fast variation of 0.27 mag in 
2.6 hr occurred in the R band; and that in any band this source 
has presented only small, short-term variations, within a total 
range of less than 0.4 mag (Viflata et al. 1997). The BVR hght 
curves of 3C 454.3 from November 1994 to December 1997, 
were presented by Villata et al. (2001). 

The VRI hght curves are shown in Fig. 11. During our ob- 
serving runs, the overall magnitude variations were AV = 0.79, 
AR - 0.93, and Al = 0.95. The maximum R magnitude 
14.75 ± 0.01 was found on October 30 (JD=2452942.662), 
while the minimum 15.67 + 0.01 was on September 29 
(JD=2452911.810). On September 29 and October 1, the 
source was stable at R = 15.67 ± 0.01 (V = 16.11 ± 0.01, 
/ = 15.15 + 0.01). The object became brighter in late October. 
It was detected at R ~ 15.1 on October 26, and it gradu- 
ally brightened over the following four days, until it reached 
R ~ 14.8 on October 30. On November 21 and 22, it was sta- 
ble atR ~ 15.0, but faded to R ~ 15.2 on November 23. In 
December, it started atR= 15.45 ± 0.01 on December 23, then 
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Fig. 6. Intraday variations of S5 0716+71 on 10 December 2003, in the B band (left panel) and in the V band (right panel). All 
data are from the 1.0 m Mt. Lemmon telescope. 
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Fig. 7. Light curves of PKS 0735+17 in the V, R, I bands. The symbols are the same as in Fig. 1 . 



brightened slightly on December 24. Finally, an R magnitude 
of 15.13 + 0.02 was reached on December 25, as recorded by 
the Sobaeksan Observatory. 



4. Spectral variability 

In this section, we investigate the relationship between the 
spectral changes and the flux variations. The plots of the colour 
indices V - I and R magnitudes for each object are given in 
Fig. 12. In most cases, the colour indices were calculated by 
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Fig. 8. Light curves of OJ 287 in the V, R, I bands. The symbols are the same as in Fig. 1. 
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Fig. 9. Light curves of 3C 345 in the V, R, I bands. The symbols are the same as in Fig. 1 . 



coupling the data taken by the same instrument within 10 min 

intervals. In Table 4, we report the principal results of our re- 
search. Column 1 gives the object name, Col. 2 the variability 
range in R magnitude, Col. 3 the number of data points used 
to calculate the colour indices V - I, Col. 4 the linear Pearson 
correlation coefficient between the V - / and R magnitudes, and 
Col. 5 the probabiUty that no correlation is present. 



From Fig. 12, it can be seen that there are no uniform trends 

of the colour index changes with the source brightnesses, and 
the significance of the dependence of the variation of the colour 
index on the source brightness (e.g. as shown by the correlation 
coefficient r) is not related with the overall variability range in 
R magnitude, nor with the number of data points. Interestingly, 
we find that the colour indices of two out of three FSRQs, PKS 
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Fig. 10. Light curves of BL Lac in the V, R, I bands. The symbols are the same as in Fig. 1. 
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Fig. 11. Light curves of 3C 454.3 in the V, R, I bands. The symbols are the same as in Fig. 1. 



0420-01, and 3C 454.3, tend to be redder when the sources are 
brighter. Especially, we find a significant anti-correlation be- 
tween the V-I colour index and R magnitude for 3C 454.3. The 
Unear Pearson correlation coefficient is found to be r = -0.837, 
with a probability of 1 .07 x 10 that no correlation is present. 
The anti-correlation in PKS 0420-01 is relatively weak, but 
still significant; however, no significant correlation is found for 



FSRQ 3C 345. These anti-correlations imply that the spectrum 

becomes redder (steeper) when the source is brighter. However, 
this result is opposite to the common colour change tendency 
in blazars, that the spectrum becomes flatter when the source 
brightens (Ghisellini et al. 1997; Ghosh et al. 2000; Clements & 
Carini 2001; Raiteri et al. 2001; ViUata et al. 2002). In contrast 
to our FSRQs, we find that aU our BL Lac objects tend to be 
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Fig. 12. The V-I colour indices versus R magnitudes for the eight blazars. The rms colour magnitude errors are indicated by the 
vertical Unes. The symbols are the same as in Fig. 3. 



bluer (flatter) when they are brighter. We find significant posi- 
tive correlations between the colour indices and R magnitudes 
for 3C 66A, S5 0716-1-71, and BL Lac, with linear Pearson cor- 
relation coefficients of r = 0.511, r = 0.438, and r = 0.657, 
respectively, at » 99.99 per cent confidence; while only weak 
correlations are found in PKS 0735H-17 and OJ 287, at very 
low confidence levels. The different trends of the colour in- 
dex variations with brightness in FSRQs and BL Lac objects 
are somewhat unusual. However, this may indicate that there 
are different physical conditions in these two populations, al- 
though all of them were selected as being red blazars in the 
present work. Nevertheless, it should be noted that our results 
are based on limited data which span only a short period (about 
three months). Certainly, more observations are needed to rein- 
vestigate the relationship between the colour index variations 
with the source brightnesses for our sample. In particular, it 
would be necessary to confirm and/or check the possible anti- 
correlation for FSRQs with more observations, and then to fur- 
ther explore the observed differences between FSRQs and BL 
Lac objects. 

As shown in Sect. 3.3, variations on different time scales 
have been found for the most extensively monitored source in 
our campaign: S5 0716-1-71. We then explored the correlation 
between the colour indices and the source brightness on dif- 
ferent time scales. It can be seen from Table 4 that there is 
a significant correlation with linear Pearson correlation coeffi- 
cient r = 0.438, using the whole data set, which spans about 
140 days. During the rapid brightening over the 4 days of 19 
- 23 November 2003 (Fig. 4), we find a significant correlation 
with Unear Pearson correlation coefficient r = 0.656, which 



suggests that the trend of flatter spectra when brighter also 
holds for variations on time scales of days. This spectral flatten- 
ing with increasing brightness was also recognized by Villata 
et al. (2000) in the 72-hour optical light curves obtained for 
S5 0716-1-71, during flie WEBT campaign of February 1999. 
However, we do not find any significant correlations between 
the colour index and brightness during the intraday variations, 
either on 23 November 2003 (Fig. 5) or on 10 December 2003 
(Fig. 6). Ghisellini et al. (1997) detected a spectral flattening 
when the flux is higher during rapid flares; however, no corre- 
lation between spectral index and brightness level was found in 
the long-term trend. A curved trajectory of a relativistic emit- 
ting blob, or very rapid electron injection and cooling pro- 
cesses, are proposed by the authors to interpret the fast vari- 
ations. Recently, using a large database, Raiteri et al. (2003) 
found that the optical colour indices are only weakly correlated 
with brightness, and that a clear spectral steepening trend was 
observed during at least one long-lasting dimming phase (JD 
~ 2451545 - 2451630). However, different spectral behaviours 
were found on shorter time scales. 

In the case of 3C 66A, the trend of flatter spectra when 
brighter is consistent with the recent results of Vagnetti et al. 
(2003), who have found a consistent trend of B - R harden- 
ing with increasing B band flux, independent of the actual flux 
value. Ghosh et al. (2000) found that the spectral index of 
the source, between V and R bands, flattened when the source 
brightness in V remained almost constant, while decreasing in 
R, between 1997 December 15 and 1998 January 2. Also, be- 
tween 1998 January 19 and 30, the source brightness decreased, 
while the V -R colour index reddened, and that this correlated 
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Table 4. The correlations between V - 1 colour index and R magnitude, for the eight blazar objects. 



Source 


AR 


TV 


r 


Probs. 




3C 66A 


0.50 


120 


0.511 


7.34x10- 


10 


PKS 0420-01 


1.48 


55 


-0.388 


2.62x10" 


-3 


S5 0716+71 


1.37 


271 


0.438 


1.11x10" 


14 


PKS 0735+17 


0.54 


38 


0.111 


4.97x10 


-1 


OJ287 


1.40 


29 


0.354 


4.99x10- 


-2 


3C 345 


0.55 


22 


-0.009 


9.66x10- 


-1 


BLLac 


1.05 


53 


0.657 


1.33x10- 


-8 


3C 454.3 


0.93 


48 


-0.837 


1.07x10- 


16 



well with the variation of the source brightness. However, from 
an analysis of a much larger data set of an extensive multiwave- 
length WEBT monitoring campaign from July 2003 through 
April 2004, Bottcher et al. (2005) found that there is a weak 
indication of a positive hardness-intensity correlation at low 
flux states with R ^ 14.0, whereas no correlation is apparent 
at higher flux levels. These authors claimed that this might be a 
consequence of the fact that the B - R hardness actually peaks 
during the rising phase of individual outbursts. 

5. Discussion 

It is notable that our FSRQs in general, and 3C 454.3 in partic- 
ular, became redder when brighter (see Fig. 12). In other words, 
the spectrum became steeper when the object was brighter, and 
flatter when fainter. This behaviour is opposite to the common 
trends for blazars, that they become bluer when they brighten 
(Ghisellini et al. 1997; Fan et al. 1998; Massaro et al. 1998; 
Fan & Lin 1999; Ghosh et al. 2000; Clements & Carini 2001; 
Raiteri et al. 2001; Villata et al. 2002). It has been found 
by different investigators that the amplitude of the variations 
is systematically larger at higher frequencies, which suggests 
that the spectrum becomes steeper when the source brightness 
decreases, and flatter when it increases (Racine 1970; Gear, 
Robson & Brown 1986; Ghisellini et al. 1997; Maesano et al. 
1997; Massaro et al. 1998). Recent investigations on spectral 
variability also show this general trend (D'Amicis et al. 2002; 
Vagnetti et al. 2003; Fiorucci, Ciprini & Tosti 2004). From our 
investigations, the variations of the colour indices of our BL 
Lac objects indeed follow this trend, e.g. 3C 66A, S5 0716+71, 
and BL Lac. This common phenomenon may be explained in 
different ways (Fiorucci et al. 2004). It may indicate the pres- 
ence of two components that contribute to the overall emis- 
sion in the optical region, one variable (with a flatter slope orvar, 
where fy oc v""), and the other stable (with Cconst > ffvar)- It is 
also possible to explain it with a one-component synchrotron 
model: the more intense the energy release, the higher the par- 
ticle's energy (Fiorucci et al. 2004). Moreover, it could also 
be explained if the luminosity increase was due to the injec- 
tion of fresh electrons, with an energy distribution harder than 
that of the previous, partially cooled ones (e.g. Kirk, Rieger 
& Mastichiadis 1998; Mastichiadis & Kirk 2002). In addition, 
it could be due to a Doppler factor variations on a spectrum 
slightly deviating from a power law, e.g. Doppler factor vari- 
ations on a 'convex' spectrum (Villata et al. 2004). It may 
also be possible that more than one mechanism are at work, as 



pointed out by Villata et al. (2004), the variabiUty observed in 
the optical curves of BL Lac can be interpreted in terms of two 
components: a 'mild-chromatic' longer-term component and a 
'strongly-chromatic' shorter-term one, which can be likely due 
to Doppler factor variations on a 'convex' spectrum and intrin- 
sic phenomena, such as particle acceleration from shock-in-jet 
events (e.g. Mastichiadis & Kirk 2002), respectively. 

However, some evidence that the amplitudes of variations 
are not systematically larger at high frequencies has been found 
on several occasions (see, for example: Malkan & Moore 1986; 
Brown et al. 1989; Massaro et al. 1998; Ghosh et al. 2000; 
Clements et al. 2003; Ramirez et al. 2004). Based on their re- 
sults, Ghosh et al. (2000) suggested that it may not be correct to 
generalize that the amplitude of the variation in blazars is sys- 
tematically larger at higher frequency. In particular, they found 
a reddening (i.e. spectral steepening) in their optical observa- 
tions of the BL Lac object PKS 0735+17. However, our anal- 
ysis of this source shows that there is a tendency of blueing 
(i.e. spectral flattening) as the brightness increases, although 
the correlation between the colour indices and brightnesses is 
rather weak (see Fig. 12 and Table 4). Moreover, Ghosh et al. 
(2000) found that the spectral slope of AO 0235+164 remained 
almost constant when its brightness increased. They concluded 
that these characteristics cannot be described simply by energy 
losses in a pure synchrotron mechanism scenario. 

The trend that FSRQ PKS 0736+017 was redder when 
brighter was noted by Clements et al. (2003). They concluded 
that the phenomenon appeared to be more related to the nature 
of the variation, than to the host galaxy or the source bright- 
ness. Recently, Ramirez et al. (2004) detected the same pe- 
culiar tendency to redden with increased brightness in PKS 
0736+017, throughout their observations. In addition, the anal- 
ysis of their data, and of the data reported by Clements et 
al. (2003), suggested two varying modes. At low flux levels, 
small changes in flux correspond to large changes in the spec- 
tral slope, while much less pronounced spectral changes corre- 
spond to a high brightness state. In both cases, the object red- 
dens when it brightens. 

FSRQs usually show strong emission lines, and a thermal 
contribution that may be comparable to the synchrotron emis- 
sion in the optical spectral region. Their optical emission is 
contaminated by thermal emission from the accretion disk and 
the surrounding regions. Of particular relevance is the presence 
of the so-called "blue bump", or "UV bump", which flattens 
the spectral slope in the optical region. Since the thermal con- 
tribution is larger in the blue region, the composite spectrum 
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would be flatter than the non-thermal component. Then, when 
the object is brightening, the non-thermal component has a 
more dominant contribution to the total flux, and the composite 
spectrum steepens. This scenario has been used to qualitatively 
explain the low flux level data of PKS 0736+017 (Ramirez et 
al. 2004). It can also be applied to explain the general trends of 
steeper spectra when brighter in our FSRQs. 

Contrary to the large majority of AGNs, which are charac- 
terized by optical spectra with prominent emission lines, BL 
Lac objects have quasi-featureless spectra. In fact, by the defi- 
nition of this class of AGNs, the line equivalent widths should 
be very small. As originally proposed by Blandford & Rees 
(1978), the weakness of the spectral lines is most probably 
due to the fact that the underlying non-thermal continuum is 
boosted by the relativistic beaming of a jet pointing in the ob- 
server's direction. In addition, the central ionizing luminos- 
ity of BL Lac objects may be relatively too weak to produce 
prominent emission lines. Consequently, the thermal contribu- 
tion is rather small in the optical spectral region, as compared 
to the synchrotron emission, even when the source is in a low 
state. In actuality, Vagnetti et al. (2003) found that BL Lac ob- 
jects and quasars are clearly segregated in the a-fi plane (see 
their Fig. 3), where a is the average spectral slope, and p is 
a spectral variability index. The lower average spectral index 
of their eight BL Lac objects indicates the absence, or lower 
relative weight, of the thermal blue bump component. The au- 
thors further proposed that the segregation in the a-fi plane is a 
consequence of the different emission mechanisms in the opti- 
cal band: synchrotron in the case of BL Lac objects, and ther- 
mal hot spots on the accretion disk in the case of QSOs. They 
found that a simple model representing the variability of a syn- 
chrotron component could explain the spectral changes of their 
BL Lac objects. Our BL Lac objects were selected to be Low- 
energy-peaked BL Lacs (LBLs); thus, we expect to observe 
them in the descending part of the spectral power distribution, 
blueward of the peak frequency. The strong variability of BL 
Lac objects is generally attributed to the synchrotron emission. 
It is thus not surprising that they generally follow the common 
trend for blazars, i.e. bluer when brighter, as we observed for 
our BL Lac objects. 

Although the thermal contribution may in general be rela- 
tively large in the optical spectral region in FSRQs, the amount 
of this contribution can be quite different from source to source. 
It is conceivable that the thermal contribution is larger in 3C 
454.3, as compared to 3C 345, in which only a weak correla- 
tion is found. It might be interesting to estimate the thermal 
contribution in these sources, to better understand the tendency 
of the colour index to vary with brightness. Nevertheless, it is 
not easy to distinguish the thermal from the non-thermal com- 
ponents using optical data alone. It is essential to obtain simul- 
taneous observations completely covering at least the near-IR- 
to-optical part of the spectrum, to better understand the spec- 
tral variability of FSRQs and/or blazars. There is no doubt that 
more observations of blazars in general, and those blazars that 
show reddening as they brighten in particular, are needed to in- 
vestigate the relationship between colour index and/or spectral 
slope and brightness in more detail. 



6. Summary 

We present the results of our monitoring of the flux variability 
of eight red blazars from 2003 September to 2004 February. 
The main results can be summarized as follows: 

1. All sources showed strong variabilities with amplitudes 
> 0.5 mag during our 3-4 month monitoring period. Variations 
with amplitudes of over 1 mag are found in four sources. In the 
extreme cases, i.e. PKS 0420-01 and OJ 287, large variations 
with amplitudes of over 1 .4 mag were observed. 

2. Intraday variations were found in S5 0716+71. 
Variations with amplitudes > 0.15 mag were detected on 
November 23 and December 10, 2003. A rapid rising rate of 
0.002 mag per min was found on December 10, 2003. We found 
a mean behaviour of a rapid linear brightening with a rising rate 
of ~ 0.2 mag per day over 4 days, during November 19-23, 
2003. 

3. We found that two out of three FSRQs tend to be redder 
when the sources are brighter; and, conversely, all our BL Lac 
objects tend to be bluer. There is a significant anti-correlation 
between the V -R colour index and R magnitude for 3C 454.3, 
which is opposite to the common colour change tendency in 
blazars. For our 3 BL Lac objects, 3C 66A, S5 0716+71, and 
BL Lac, there are strong positive correlations, while only very 
weak coiTelations were found for the remaining 2 BL Lac ob- 
jects, PKS 0735+17, and OJ 287. 

4. We propose that the different relative contributions of the 
thermal versus non-thermal radiation to the optical emission 
may be responsible for the different trends of the colour index 
with brightness in FSRQs and BL Lac objects. 
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